Cosmological Fluctuations from Infra-Red Cascading During Inflation 
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We propose a qualitatively new mechanism for generating cosmological fluctuations from inflation. 
The non-equilibrium excitation of interacting scalar fields often evolves into infra-red (IR) and 
ultra-violet (UV) cascading, resulting in an intermediate scaling regime. We observe elements of 
this phenomenon in a simple model with inflaton (j) and iso-inflaton x fields interacting during 
inflation via the coupling g^{(f> — (po)^x'^- Iso-inflaton particles are created during inflation when 
they become instantaneously massless at = 0o, with occupation numbers not exceeding unity. 
Previous studies have focused on the momentary slowing down of the condensate by back- 
reaction effects. Here, we point out that very quickly the produced x particles become heavy 
and their multiple re-scatterings off the homogeneous condensate <j){'t) generates Bremsstrahlung 
radiation of light inflaton IR fluctuations with high occupation numbers. The subsequent evolution 
of these IR fluctuations is qualitatively similar to that of the usual inflationary fluctuations, but 
their initial amplitude is different. The IR cascading generates a bump-shaped contribution to 
the cosmological curvature fluctuations, which can even dominate over the usual fluctuations for 
> 0.06. The IR cascading curvature fluctuations are significantly non-gaussian and the strength 
and location of the bump are model-dependent, through g^ and (f)o- The effect from IR cascading 
fluctuations is significantly larger than that from the momentary slowing-down of 4>{t). With a 
sequence of such bursts of particle production, the superposition of the bumps can lead to a new 
broad band non-gaussian component of cosmological fluctuations added to the usual fluctuations. 
Such a sequence of particle creation events can, but need not, lead to trapped inflation. 
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I. INTRODUCTION 

In addition to the standard mechanism for generat- 
ing cosmological perturbations during inflation from the 
vacuum fluctuations of the inflaton field [l| , there are al- 
ternative mechanisms including modulated fluctuations 
(inhomogeneous preheating) 0,0] and the curvaton 4], 
both of which are based on the vacuum fluctuations of 
iso-inflaton fields during inflation. In this paper we pro- 
pose a new and qualitatively different mechanism for gen- 
erating cosmological fluctuations during inflation. 

Physical processes during inflation may leave their im- 
print as features in the cosmological fluctuations. These 
can, in principle, be observed if they fall in the range of 
the wavelengths between 10^ Mpc and 100 Kpc, which 
corresponds to about ten e-folds during inflation. There 
may also be signatures such as additional features at the 
horizon scale or potential observables on much smaller 
scales. Relevant dynamical models were studied in the 
early days of the inflationary theory, for example the 
model with phase transitions during inflation yielding as- 
sociated features in the cosmological fluctuations 5, 6, 7] 
(see also [8] ) . There the time-dependent dynamics of the 
inflaton field (f> can trigger a phase transition in the iso- 
inflaton X field. The growth of x inhomogeneities induces 
curvature fiuctuations on scales leaving the horizon at the 
moment of the phase transition. 

Recently, several studies d, H, [l^, [HI [12] considered 
features in the cosmological fluctuations from the effect of 
particle creation during inflation, which can be modeled 



by the simple interaction 
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with some value of the scalar fleld 0o which the rolling 
(l){t) crosses during inflation that must be tuned to give 
a signal in the observable range of e-folds. There are dif- 
ferent motivations for the model The early study 
[l3l | introduced the possibility of slowing down the fast 
rolling inflaton using particle creation via the interac- 
tion ll]). Imagine there are a number of field points 4>oi, 
i — 1,2, ■ ■ ■ ,n, where the iso-inflaton field becomes mass- 
less and X particles are created. The produced x particles 
are diluted by the expansion of the universe, however, the 
back-reaction effect from multiple bursts of particle cre- 
ation may slow down the motion of (j) sufficiently to allow 
for slow-roll inflation. This is called trapped inflation. A 
more concrete string theory realization of trapped infla- 
tion, based on the sequence of I?3 branes interactions, 
was discussed in [l^]. The work [l2], which is compli- 
mentary to this study, provides a detailed realization 
of trapped inflation in the context of the string theory 
model (13 • 

The instant of x-particle creation and the slow-down of 
the rolling inflaton shall generate a feature in the power 
spectrum of scalar curvature fluctuations from inflation 
P(;{k). This was noticed in |9|, where the features in the 
power spectrum were estimated from the simple-minded 

formula P({k) ^ ('^) • inflaton slow-down was 
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described by the mean-field equation 



(2) 



The vacuum expectation value (VEV) (x^) can be cal- 
culated with the analytic machinery of particle creation 
with the coupling ([1]) , which was developed in the theory 
of preheating after inflation The QFT of X parti- 

cles interacting with the time-depended condensate (j)(t) 

deals with the eigenmodcs Xkit)e^'''^, where the time- 
dependent mode function obeys an oscillator-like equa- 
tion in an expanding universe 



Xfe - , (3) 



with time-dependent frequency ujk (t) . When (j){t) crosses 
the value (f>o, the Xk mode becomes massless and ujk{t) 
varies non-adiabatically. Around this point {(j){t) — (j)Q) w 
00 (t — to), where t = tn is corresponding time instant. 
With this very accurate [l6| approximation, one can solve 
the equation ^ analytically to obtain the occupation 
number of created x particles 



rik = exp 



"fcT 



(4) 



presuming that k^, > H. The latter condition requires 
coupling constant g > iJ^/|0o| ^ 10^'*. It is useful to 
note that, independent of the details of V{(p) and (p{t), 
the scale fc* can be related to the naively estimated am- 
plitude of vacuum fluctuations as = ^ g / {2ttV^^'^). 

Thus K/H - 30 if Vl'"^ = 5 x 10"^ as suggested by the 
CMB and the coupling is ~ 0.1.^ 

The VEV (x^): which controls the back-reaction on the 
homogeneous field (j){t), can be calculated from (|4]) and 



estimated as {x^] 
for 6 > 



_ r d-'k 
— J TTTv^ 



(2ir)^ iXfcl ~ J (27r)^ Uk ~ g\ip- 

tiQ. Substitution of this results back into ([2|) 
gives expected velocity dip of (f){t) and, correspondingly, 
a bump in the power spectrum P(;{k). In Fig. [T] we illus- 
trate this velocity dip for the model ^ with g^ = 0.1. 

The calculation of curvature fluctuations in the model 
P]) was re-considered in where the linearized equa- 
tions of motion for the quantum fluctuations Scf) cou- 
pled with the metric fluctuations were treated again in 
the mean-field approximation, using (x^) to quantify the 
back-reaction. This study shows that the bump in the 
curvature power spectrum is the most prominent part of 
an otherwise wiggling pattern. Similar to us, the work 
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FIG. 1: \(f>\/{Mpm) plotted against mt for g = 0.1 (where 
m = 1^0,^ is the effective inflaton mass). Time t = corre- 
sponds to the moment when cj> = cj>Q and x-particles are pro- 
duced copiously. The solid red line is the lattice field theory 
result taking into account the full dynamics of re-scattering 
and IR cascading while the dashed blue line is the result of a 
mean field theory treatment which ignores re-scattering pdl ]. 
The dot-dashed black line is the inflationary trajectory in the 
absence of particle creation. 



[12| | further refined the calculation of the curvature per- 
turbation in this model, going beyond the mean-field 
treatment of 6? 




X 



FIG. 2; Re-scattering diagram. 

In a parallel development, scalar fields interactions of 
the type ([1]) are the subject of studies in non-equilibrium 
QFT and its application to the theory of preheating af- 
ter inflation, as we mentioned above. Although we study 
particle production during inflation (as opposed to dur- 
ing preheating, after inflation) there are many similari- 
ties. For example, in the case of parametric resonant pre- 



^ We are assuming that supersymmetry protects the inflaton po- 
tential from radiative corrections at t = to- An explicit realiza- 
tion of the type of coupling we are interested in, based on global 
Af = 1 supersymmetry, has been provided in IITII . For string 
theory models the reader is referred to llCt [T^ . If4ll . 



^ Below we use QFT methods to study correlators of inhomoge- 
neous fluctuations Sif) induced by inhomogeneities. Ref. [T^ 
considers the effect induced by quantum mechanical fluctuations 
of the total particle number n^- Owing to the relationship be- 
tween and Tij^., our calculations below capture this effect. 
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heating due to the oscillating inflaton background (j){t), 
X particles are created in successive bursts whenever (/){t) 
crosses zero and the % particles become instantaneously 
massless. This leads to huge occupation numbers of the 
created x fluctuations. On the other hand, in the sce- 
nario described above there is only a single burst of par- 
ticle production and the resulting x occupation number 
(HI is always less than unity. 

The full dynamics of interacting scalars during pre- 
heating also includes not only bursts of particle pro- 
duction but also re-scattering effects where d(f) fluctua- 
tions (particles) are created very quickly due to the in- 
teraction of created x particles with the condensate (j){t) 
[TgI [25I [26| . The diagram for this process is illustrated 
in Fig. O The Scj) particles produced by re-scattering are 
far from equilibrium and evolve towards an intermedi- 
ate regime which is well described by the scaling "turbu- 
lent" solution. To understand the dynamics, one needs 
to use lattice numerical simulations of time-evolution of 
the classical scalar fields based on the LATTICEASY ^ 
or DEFROST fl^ codes, designed for this purpose. The 
turbulent regime of interacting scalars was investigated 
in several recent works. The papers [13, |2l| used nu- 
merical simulations to demonstrate the scaling regime in 
the model of self-interacting classical scalar Xcj)'^. The 
papers [2^ [2^ show numerically the scaling solution for 
the fully QFT treatment of the same model, and advo- 
cate the new regime, the non-thermal fixed point, which 
may be asymptotically long (before the system evolves, 
if ever, to another fixed point: thermal equilibrium). 

In this paper we will study in detail the back-reaction 
of X particles produced during inflation on the inflaton 
fleld, resulting in Bremsstrahlung radiation and IR cas- 
cading of 6(p fluctuations. Our results will also apply to 
the early stages of rescattering in preheating after infla- 
tion (this is so because the time scale for re-scattering is 
short and hence the results are insensitive to the expan- 
sion of the universe). 



II. NUMERICAL STUDY OF RE-SCATTERING 

To study the creation of Scj) fluctuations by re- 
scattering of produced x particles off the condensate <j){t) 
in the model ([T]) we have adapted the numerical DE- 
FROST code for the problem of a single burst of instan- 
taneous particle creation during inflation.'^ To run the 
classical scalar field simulation, we must first choose the 
appropriate initial conditions. The field x on the lat- 



Since the production of long wavelength (50 modes is so energeti- 
cally inexpensive, a major requirement for successfully capturing 
this effect on the lattice is respecting energy conservation to very 
high accuracy. In our modified version of DEFROST energy con- 
servation is respected with an accuracy of order 10~*, compared 
to 10"^ - 10-5 obtained using previous codes. A minimum ac- 
curacy of roughly lO"'' is required for this problem. 



tice is modeled by the random gaussian field realized as 
the superposition of planar waves Xk{t)e^''^ with random 
phases. The initial conditions for the models Xk{t) are 
chosen to emulate the exact quantum mode functions cor- 
responding to the physical occupation number ^ (see 
appendix A for more details) while ensuring that the 
source term for the Scj) fiuctuations turns on smoothly at 
t — 0. The box size of our 512'^ simulations corresponds 
to a comoving scale which initially is |^ ~ 3 times the 
horizon size 1/H, while k^, = 60y/gH. We run our sim- 
ulations for roughly 3 e-folding from the initial moment 
to when x-particles are produced, although a single e- 
folding would have been sufficient to capture the effect. 
We are interested in the power spectrum of inflaton fluc- 
tuations Ptf, = fc'^|(50p/(27r^), and also the number den- 
sity of inflaton fluctuations n^{k) = ^ (^^^r^ + I'^^fep) 

(where wc introduce the notation Qk = \Jy.tt»t> + for 
the inflaton frequency). For the sake of illustration we 
have chosen the standard chaotic inflationary potential 
y = m^cjP-j^ with = IQ-^Mp and 0o = 3.2Mp, how- 
ever, our qualitative results will be independent of the 
choice of background inflation model** and, in particu- 
lar, are applicable to trapped inflation. We have con- 
sidered three different values of the coupling constant, 
5^ = 0.01,0.1,1, although we focus most of our atten- 
tion on the case — 0.1. Fig. [3] shows time evolution of 
the re-scattered inflaton power spectrum P,f,{k) for three 
different time steps, while Fig. 3] shows the correspond- 
ing evolution of the particle number density n^(k). In 
Fig. [5] we illustrate the dependence of our results on the 
coupling constant g^. 

In Fig. [3] we see clearly how multiple re-scatterings lead 
to a cascading of power into the IR. These re-scattered 
inflaton perturbations are complimentary to the usual 
long-wavelength inflaton modes produced by quantum 
fluctuations. As long as g^ > 0.06 the re-scattered power 
spectrum outside the horizon comes to dominate over the 
usual vacuum fluctuations within a single e-folding. At 
much later times the IR portion of the power spectrum 
remains frozen while the UV portion is damped out by 
the Hubble expansion. The effect of IR cascading on the 
power spectrum is much more significant than features 
which are produced by the momentary slowing-down of 



* The choice of background inflationary potential will alter the 
functional form of <ji{t), however, all the dynamics of re-scattering 
occurs within a single e-folding of the moment when (p = (j>o. 
Hence, for any inflation model it will be a good approximation 
to expand (j>{t) = <po + vt and our results should depend only on 
the velocity v = <j>{t = 0) (assuming this to be non-zero) which 
is related to the Hubble scale and the observed amplitude of 
the curvature perturbation. The claim of model independence is 
borne out by explicit analytical calculations in the next section. 
There we find that the dominant contribution to P/^ from IR 
cascading depends sensitively only on the ratio fc* /H which, as 
we have shown, is determined entirely by the coupling g^. 
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FIG. 3: The power spectrum of inflaton modes induced by re- 
scattering (normalized to the usual vacuum fluctuations) as 
a function of ln(fc/fc*), plotted for three representative time 
steps in the evolution, showing the cascading of power into 
the IR. For each time step we plot the analytical result (the 
solid line) and the data points obtained using lattice held 
theory simulations (diamonds). The time steps correspond to 
the following values of the scale factor: a = 1.03, 1.04, 2.20 
(where a = 1 at the moment when = ^o). By this time the 
amplitude of fluctuations is saturated due to the expansion of 
the universe. The vertical lines show the range of scales from 
our lattice simulation. 




FIG. 4: Physical occupation number Uk as a function of 
\n{k/K) for g'^ = 0.1. The three curves correspond to the 
same series of time steps used in Fig. |3] and demonstrate the 
growing number of long wavelength inflaton modes which are 
produced as a result of IR cascading. Because the same x- 
particle can undergo many re-scatterings off the background 
condensate <?i(t), the d(j) occupation number is larger than 
the initial x particle number (for =0.1 one can achieve 
n^{k) ~ 30 even though initially n^(k) < 1). When g^ — 0.06 
the IR S(j) occupation number exceeds unity within a single e- 
folding. The yellow envelope line shows the early onset of scal- 
ing behaviour associated with the scaling turbulent regime. 



the background (t>{t).^ 

Fig. [T] illustrates the impact of re-scattering on the dy- 
namics of the velocity of the background field. The evo- 
lution of (/)(t) including re-scattering is not changed sig- 
nificantly (as compared to the mean field theory result), 
which show the energetic cheapness of IR cascading. 

The long- wavelength inflaton fluctuations produced by 
IR cascading are non-gaussian. This is illustrated in 
Fig. [6] where we study the probability density function 
and compare to a gaussian fit. 



To avoid confusion: here we use "cascading" to refer to the dy- 
namical process of building up S<j> fluctuations in the IR. If the 
universe were not expanding, a scaling turbulent regime would be 
established. Here we see this scaling regime only in an embryonic 
form, see the envelope in Fig. |4] 



III. ANALYTICAL THEORY OF 
RE-SCATTERING 

We now develop an analytical theory of this effect. 
Here we provide only a cursory discussion, the reader 
is referred to Appendices A and B for a detailed exposi- 
tion and technical details of the calculation. At leading 
order the physics of re-scattering (see Fig. [J) is described 
by the equation 

(50 + iH5^ - \^^5(t) + m^Scf) = -g^ [0(<) - M ^ , (5) 

where we introduce the notation w? — Vc/icii for the ef- 
fective inflaton mass (hence we are not assuming a back- 
ground potential of the form to^(/)^/2 in this section, only 
that ^ in the vicinity of the point (j) — (f)o)- The 
solution of ^ consists of two components: the solution 
of the homogeneous equation which simply corresponds 
to the usual vacuum fluctuations produced during in- 
flation and the particular solution which is due to the 
source term. We will focus our attention on this latter 
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FIG. 5: The dependence of the power spectrum on the 
coupling g^. The three curves correspond to for = 
0.01, 0.1, 1, evaluated at a fixed value of the scale factor, 
a — 2.20. We see that even for small values of g^ the in- 
flaton modes induced by re-scattering constitute a significant 
fraction of the usual vacuum fluctuations after only a single 
e-folding. 



solution which, physically, corresponds to re-scattered in- 
flaton perturbations. Since the process of IR cascading 
takes less than a single e-folding, we can safely neglect the 
expansion of the universe when studying analytically the 
particular solution of ([5]) . (In all of our lattice simulations 
the inflationary expansion of the universe is taken into 
account consistently.) Solving for the particular solution 
(5(/ifc of ^ and defining the re-scattered power spectrum 
in terms of the QFT correlation function in the usual 
manner we arrive at an expression for in terms of the 
c-number mode functions Xk which obey equation Q.^ 
The result is 



J dt'dt"t't" sin [^kit - t')] sin \nk{t - t")] 



X / d^k'x,.,.{t')xl.k,{t")xk'{t')xl,{t") (6) 



for the (50-particle 



where again we have fife 
frequency. 

To evaluate this power spectrum we need an expression 
for the solutions of ^ in the regime of interest. Let us 
choose the origin of time so that i = corresponds to 
the moment when <j> = (j>o. At the moment t = the 
parameter \ujk\/uj'^ is order unity or larger and ujk varies 
non-adiabatically. At this point Xfe modes are produced 
in the momentum band k ^ ki,. However, within a time 
At ~ fc~^ (which is tiny compared to the Hubble time 
H^^) the X particles become extremely heavy and their 
frequency again varies adiabatically. At times t ^ k~^ 
we can safely approximate Wfe = yjk"^ + k^t^ = k^t for 
the modes of interest and Xk takes the simple form 




g-i(fe,t)V2 ^+iik,tf/2 

Xk{t) = Vl + nk — ; 7^ i\/rik 



k^V2i 



kW2t 



, (7) 



10'5(|)/M^ 



where the occupation number was defined in The 
factors VI + i^k, —iy/nk are the Bogoliubov cocfflcicnts 
while the factors proportional to e^'*^'^**-' come from 
the positive and negative frequency adiabatic mode so- 
lutions As we see, very quickly after t = 
the X particles become very massive and their mul- 
tiple re-scatterings off the condensate (j){t) generates 
Bremsstrahlung radiation of IR S(j) particles. 

We have computed the full renormalized power spec- 
trum analytically in closed form and the result is pre- 
sented in equation (jA-17[) . This formula is used for all 
of our figures. Since the exact analytical result is quite 
cumbersome, it is useful to consider the following repre- 



FIG. 6: Probability density function of S4> for = 1 at a 
series of different values of the scale factor, a. The dashed 
curve shows a Gaussian fit at late time a = 6.2. 



® We are only interested in connected contributions to the correla- 
tion functions, which is equivalent to subtracting the expectation 
value from the source term in JSjl: ^ ~ ix^)- Thus, our 
re-scattered inflaton modes are only sourced by the variation of 
X^ from the mean {x^)- 
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sentative contribution to ([6]): 



9' 



k^k 



3 r 



32V2 7r5 



1 — COs(f2fet) 



(8) 



which captures the properties of the full analytical solu- 
tion. In particular, the simple expression ([5]) nicely de- 
scribes the IR cascade. The spectrum has a peak which 
initially (near t ^ k~^) \s close to ki,. As time evolves the 
peak moves to smaller-and-smaller k as power builds up 
in the IR. From ([5]) we see that modes with fifct < 1 gain 
power as Pif>{k) ~ t^. For a given fc-mode the growth of 
the power spectrum saturates when Vlkt ~ 1, however, 
the cascade still continues at some lower k. If we had 
m = then the cascade would continue forever, other- 
wise formula ([8]) predicts that the growth of the spectrum 
saturates at t ~ when the peak has reached k ^ m. 
After this point the character of the IR cascade is ex- 
pected to change, however, our analytic calculation is no 
longer reliable because m H and we have neglected the 
expansion of the universe. Notice from equation ([5]) that 
the value of /H'^ at the peak (which is the main ob- 
servable signal) is fixed by the ratio k^,/ H , which in turn 
depends only on the coupling (there is some depen- 
dence also on = V^^^, however, this is subdominant 
since ^ for any inflation model). This observa- 
tion confirms our previous claim that the dynamics of IR 
cascading are largely insensitive to the choice of back- 
ground inflation model. 



IV. DISCUSSION OF CURVATURE 
FLUCTUATIONS FROM IR CASCADING 

Any inflaton fluctuations 5(j), independently on their 
origin, evolve qualitatively similarly during inflation. 
When their physical wavelength is smaller than the Hub- 
ble radius Sip is oscillating while their amplitude 
is diluted as 1/a. As far as the wavelength exceeds the 
Hubble radius, the amplitude of Sip freezes out. Fluc- 
tuations of 54) induce the curvature metric fluctuations. 
Inflationary expansion of the universe further stretch the 
wavelengths of the fluctuations frozen outside the hori- 
zon, making them potentially of the cosmological scales, 
depending on the wavelength. The inflaton fluctuations 
produced by the IR cascading, therefore, are the poten- 
tial sources for observable curvature fluctuations. To cal- 
culate curvature fluctuations generated by the IR cascad- 
ing, we have to solve self-consistent system of linearized 
Einstein equations for metric and the fields fluctuations. 
For example, the (0, 0) linearized Einstein equation for 
our model reads 



<5G° = |J(5r„"(0)+<5T°(x)) 



(9) 



where (5Gq is the perturbed Einstein tensor and in the 
r.h.s. 6TQ{(t)) corresponds to the fluctuations of the infla- 
ton energy density, containing familiar terms linear with 




100 



50 



= 0.1 




-2 -1 
In (k/k*) 



FIG. 7: The top panel shows a comparison of curvature fluc- 
tuations from different efi'ects. We see the dominance of fluc- 
tuations produced by IR cascading over the wiggles induced 
by the momentary slowing-down of the inflaton. For illustra- 
tion we have taken — 0.1, but the dominance is generic for 
all values of the coupling. The red solid line is the IR cascad- 
ing curvature power spectrum, while the blue dashed line is 
the result of a mean field treatment. (The vertical lines show 
aH at the beginning of particle production and after ~ 3 
e-foldings.) The bottom panel shows the curvature power 
spectrum resulting from multiple bursts of particle produc- 
tion and IR cascading. Superposing a large number of these 
bumps produces a broad-band spectrum. 



respect to S(f>, like (pScj) etc. Second term corresponds to 
the contribution from x particles 



M'x'- < T^ix) > ■ 

(10) 

Although this expression is bi-linear w.r.t. x, it turns 
out taking STq{x) into account is important. To be- 
gin the investigation the equation Q, it is convenient 
to use its Fourier transformation. While the Fourier 
components of JTq (0) contains linear terms of 6(pk, the 
Fourier transform of 6Tq(x) contains convolutions like 
~ / ifik'xpxt p. etc. As a result, despite the fact 
that X particles amplitude is peaked at fc '-^ fc^, this type 
of convolution gives significant contribution at small k, 
which are of interest for the theory of generation of cos- 
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mological fluctuations. Preliminary estimations based on 
the analytical formulas for Xk involved in the convolution, 
show that contribution of 6Tq{x) is at least the same or- 
der of magnitude as STq{(P). Rigorous treatment of the 
curvature fluctuations in our model is therefore rather 
complicated, and will be leave it for separated project. 

Since both terms in r.h.s. in ([9]) are the same order of 
magnitude, here for the crude estimations we will use the 

simple-minded formula ~ ('^) ' The curvature fluc- 
tuations generated by the IR cascading, are illustrated 
in the top panel of Fig. [T] The curvature fluctuations 
from the instance of the IR cascading has the bump- 
like shape within the interval of the wavelength, roughly 
corresponding to one e-folding. They are significantly, 
by orders of magnitude, dominated over the fluctuations 
generated by the momentary slowing down of (j){t) . If we 
pick up the background inflationary model to be chaotic 
inflation with the standard quadratic potential, the ratio 
of the power spectra from IR cascading and the standard 
fluctuations is estimates as Pjr/Ps ^ 700 x g'^-^. Thus, 
depending on the coupling 5^, the IR bump can domi- 
nate (for g"^ > 0.06) over the standard fluctuations, or 
just contribute to them for smaller g^. 

Suppose that we have a sequence of the particles cre- 
ation events at different moments toi, « = 1, 2, 3, .... Each 
of those events generate, through IR cascading, corre- 
sponding bumps in the spectrum, as illustrated in the 
lower panel of Fig. [T] Depending on the density of ioi mo- 
ments, superposition of such IR bumps results in broad 
band contribution to the curvature power spectrum. 

We also estimated non-gaussianity of Sip fluctuations 
from IR cascading. They are quite signiflcant, we esti- 
mate the non-gaussianity parameter /jv^ ~ 2 x 10* g'^ '^^. 
Therefore the non-gaussian signal from individual bump 
can be strongly non-gaussian. In the model with mul- 
tiple instances of particle creations, the broad-band IR 
cascading fluctuations dominated over the standard fluc- 
tuations, apparently, are ruled because of the strong non- 
gaussianity. However, the broad-band IR cascading fluc- 
tuations can be considered as additional subdominant 
component to the standard fluctuations. In this case the 
non-gaussianity of the net curvature fluctuations can be 
acceptable but different from that of the standard fluctu- 
ations alone. We leave a detailed discussion of the non- 
gaussianities produced during IR cascading (and their 
observability) to future studies. Notice that this type of 
nongaussianity, which is localized over a narrow region 
of scales, is similar to what was obtained in [27|. 

Another important parameter of the IR cascading fluc- 
tuations is the wavelength of the bump, which depend 
on the value of (/>o. there are interesting possibilities to 
consider them at small CMB angular scales (small-scale 
non-gaussianity?), at scales of galaxies, at near the hori- 
zon scales (CMB anomalies at large scales?). We leave 
all of these possibilities for future discussion. 



V. SUMMARY AND CONCLUSIONS 

We find the following new results for interacting scalars 
during inflation in the model ([1]). 

i) In the early stages of re-scattering, when the back- 
reaction can be treated linearly, the spectrum of infla- 
ton fluctuations 6(p and the corresponding particle num- 
ber density n^{k) can be rigorously calculated with QFT 
with the diagram in Fig.[2l We perform such QFT calcu- 
lations and compare with lattice simulations of the clas- 
sical field dynamics. The results are highly compatible 
with each other, even well into the late time nonlinear 
regime. This signals the dominance of Fig. [2] in the dy- 
namics of re-scattering, while the analytic estimate gives 
a handy fitting formula. 

ii) While the stationary scaling turbulent solution for the 
scalar fields after preheating was established previously, 
the way this regime appears dynamically was not traced 
out in detail. For our example for the first time we show 
explicitly how this scaling regime is seeded. In the ab- 
sence of expansion of the universe, this embryonic scaling 
behaviour will develop into the full turbulent regime that 
has previously been observed, however, for our purposes 
only the early stages are relevant because fluctuations 
freeze outside of the horizon. 

iii) The most unexpected result, which is of interest 
outside of the inflationary theory, is that even an in- 
significant amount of out-of-equilibrium particles with 
"n-xik) < 1, being re-scattered off the scalar field conden- 
sate, can generate IR cascade of the inhomogeneous con- 
densate fluctuations with large occupation number n^{k) 
in the IR region. This is explained by fact that multiple 
production of the IR modes is energetically cheap. 

iv) IR fluctuations of the light flelds have special sig- 
nificance in the context of inflationary theory. These 
fluctuations evolve in time similar to the evolution of 
the usual inflationary fluctuations. Their amplitude is 
oscillating while their wavelengths is inside the Hubble 
radius and is frozen out once their wavelengths exceed 
the Hubble radius H^^ . However, the amplitude of IR 
cascading fluctuations is different from that of the usual 
quantum fluctuations. Frozen fluctuations (50, regardless 
of their origin, will induce cosmological curvature fluc- 
tuations. Thus, we get a new mechanism for generating 
frozen long wavelength 5(j} fluctuations from IR cascad- 
ing. Therefore, IR cascading will lead to observable fea- 
tures in the CMB power spectrum. For generic choices 
of parameters, these re-scattered fluctuations are much 
more significant than the features induced by the mo- 
mentary slowing-down of the background (/'(t), see the 
upper panel of Fig. [T] 

v) Since the solution 5(j} of ([5]) depends nonlinearly on the 
gaussian field x, the curvature fluctuations induced by IR 
cascading will be non-gaussian. This nongaussianity is il- 
lustrated in Fig. [6l We estimate this nongaussianity to be 
significant. However, this non-gaussian signal is related 
to the IR cascading bump of the spectrum and peaks on 
the range of scales corresponding to roughly one e-folding 
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after t = to- This type of nongaussianity, which is large 
only over a small range of scales, is not well constrained 
by observation. 

vi) The strength and location of our effect is model- 
dependent (through and (po), however, the very fact 
that subtle QFT effects of interaction during inflation 
may lead to an observable effect is intriguing. 

vii) In our analysis we have focused on a single burst of 
instantaneous particle production during inflation. This 
scenario is interesting in its own right, however, our re- 
sults could also be extended in a straightforward manner 
to study trapped inflation models where there are nu- 
merous bursts of particle production; see for more 
detailed discussion. 

viii) Suppose we have a sequence of points (jjoi (* — 
1, • • • , N) where particles Xi become massless. In this 
case the curvature fluctuation profiles generated from in- 
dividual bursts of IR cascading can superpose to form 
a smooth spectrum of cosmological fluctuations, see the 
lower panel of Fig. [T] This provides us with a new mech- 
anism for generating long wavelength curvature fluctua- 
tions during inflation from IR cascading. The amplitude 
and non-gaussianity of these curvature fluctuations will 
depend on the coupling, g'^. These fluctuations are in- 
teresting on their own, although they may generate too 
much non-gaussianity. They also can be considered as an 
extra component of the standard vacuum fluctuations, 
introducing an interesting non-gaussian signal to the net 
fluctuations. 

ix) The transfer of energy into fluctuations via succes- 
sive bursts of particle production can lead to trapped 
inflation. Our new mechanism of generating cosmologi- 
cal fluctuations from IR cascading can, but need not, be 
associated with trapped inflation. 

x) Varying the location, strength and non-gaussianity of 
the IR cascading bump, it will be interesting to consider 
other potential implication to the cosmological fluctua- 
tions, e.g. their impact on the horizon scale fluctuations 
or on small scale fluctuations where they might effect pri- 
mordial Black Hole formation or the generation of grav- 
itational waves. 



APPENDIX A: Analytical Theory of Re-scattering 

In this appendix we develop an analytical theory of 
re-scattering which is in good agreement with the re- 
sult of fully nonlinear lattice field theory simulations. 
As usual we split the infiaton field into a classical ho- 
mogeneous component and quantum inhomogeneities as 
(f>{t, x) = (f>{t)+S(j){t, x) such that {(f>{t, x)) = and we 
further suppose that (x(i,x)) = 0. Since IR cascading 
occurs within a single e-folding we can safely neglect the 
expansion of the universe. However, there is no obstruc- 
tion to consistently including this effect [1^. 

At leading order the physics of re-scattering is de- 
scribed by equation ([5]), which corresponds to the dia- 
gram in Fig. [2] There is a correction to ([5]) corresponding 
to a diagram where two S4> particles interact with two x 
particles, however, this effect is sub-leading [l^. It is un- 
derstood that one must subtract from ([5]) the expectation 
value of the right-hand-side in order to consistently de- 
fine the quantum operators Scj) such that (Scj)) = {x) — ^■ 
Subtracting off this expectation value is equivalent to 
only considering connected diagrams when we compute 
correlation functions. 



A.l. Production of x-Particles 

To solve equation ([5]) we first require explicit expres- 
sions for the background field (j){t) and the wavefunction 
x(t,x). Let us choose the origin of time so that (j) — (po 
at t = 0. Near the moment of particle production we 
can expand (f)(t) — </)g = (pQt. The interaction term in ([T]) 
induces induces a mass for the x-field 



ml = g' m) - 0o]' = f/4>'t' = ktt^ 



(A-1) 



which vanishes at i = 0. At this moment particles will 
be copiously produced by quantum effects. 

The mode functions Xk{t) obey the following equation 



Xk{t)+iol{t)xk{t)^Q 



(A-2) 



Finally, let us return to the old story of how cosmolog- 
ical fiuctuations are affected by phase transition during 
inflation, which we discussed at the beginning of this pa- 
per. We project that our results concerning re-scattering 
and IR cascading will radically change the conventional 
picture. 



where the time-dependent frequency is 



iOk{t) = ^k^ + ml = v/fc2 + fc2(M)2 (A-3) 



The theory of equation (|A-2[) is well-studied in the lit- 
erature [lo, [iBl- As long as the frequency (jA-3[) varies 
adiabatically Icj^l/Lj^ ^ 1 the modes of x will not be 
excited and are well described by the adiabatic solution 
Xk{t) — fk(t) where we have defined 
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fk{t) 



I 



exp 



I / dt'ujk{t') 



(A-4) 



However, very close to i = 0, roughly within the interval 
—k~^ < t < +k~^, the parameter jw^j/w^ can become 
order unity or larger for low momenta k ^ k^, and Xk 
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modes within this band will be produced. The general 
solution of (|A-2[) can be written in terms of the adia- 
batic modes (jA-4| and the time-dependent Bogoliubov 
coefficients as 



Xk{t) ^ ak{t).fk(t) + Pk{t).fk(tr 



(A-5) 



where the Bogoliubov coefficients obey a set of cou- 
pled ordinary differential equations with initial condi- 
tions |afe(0")| = 1, Pk{0-) = 0. Near t = the adia- 
baticity condition is violated and (3^ grows rapidly away 
from zero as a step-like function. Very shortly after this 
burst of particle production the frequency again varies 
adiabatically and ak, (3k become constant, taking the 
following values p^ : 



akit>0) 
Pk{t>0) 



/Uk e 



nk 



(A-6) 
(A-7) 



where the physical occupation number is defined by 
The phase 5k has been computed analytically in |15j and 
depends nontrivially on k. However, since most of the 
particle production occurs for momenta fc <, fc* it is an 
excellent approximation to use the simple result e*'^'= = 
— i. (We have verified that changing the relative phase 
will at most alter factors order unity in the final results.) 

We are now in a position to write out the solution for 
the Xk modes in the outgoing adiabatic regime t XL k~^. 
Since most of our interest is in IR modes with fc ;< fc* it 
is a good approximation to expand the frequency (|A-3[) 
as uJk{t) = ki,{ki,t). Using the equations (|A-6P and (jA-7P 
we can write the solution (|A-5P in the region of interest 
as 

g-j(fe*t)V2 g+i(fc,t)V2 

Xk{t) = \/l + Uk — ; — j= i^/n^ — — j^- (A-8) 



k^VTt 



A. 2. Equations for Re-scattering 

Having reviewed the solutions for (j){t) and Xk{t) we 
now turn our attention to solving Let us first briefly 
discuss our conventions for fourier transforms and mode 
functions. We write the q- number valued fourier trans- 
form of X a-s 



X(i,x) 



(Pk' 

(27r)3/2' 



(A-9) 



Because x is gaussian we can expand into c-number 
mode functions Xk (discussed above) and annihila- 



tion/creation operators Ok, a]^ as 

&) ^ a^,Xk{t) + al^xlit) 



(A-10) 



In the theory of preheating/moduli trapping without 
re-scattering the distinction between q-number fourier 
transform and c-number mode functions is not important 
because both obey the same equation of motion (equation 



(IA-2[) in the case at hand). However, once re-scattering 
is taken into account this distinction is crucial. To see 
why, note that the solution 54) of equation lU will not be 
gaussian and hence will not admit an expansion of the 
form (|H0|. 

Finally, we return to the equation for re-scattering, 
eqn. ([5|) . We can solve for the q- number fourier transform 
of 5(j) (defined analogously to (|A-9j) ) using the retarded 
Green function 



1 



(27r)3/2 JO 



dt't' sin [Vlk{t~t')] 
X I d'k'(,l_^,[t')il,[t') (A-11) 



where we have introduced the notations Vlk = Vk"^ + nfi 
for the (5(/)-particle frequency and = V,tp^ for the effec- 
tive 54> mass. Carefully carrying out the Wick contrac- 
tions yields 



X j dt'dt"t't" sin [^kSt ~ t')]s\n [Vlk^lt-t")] 
X / d^k'xk,-k'{t')xt^k'{t'')XM'{t')xUt'') (A-12) 



where the x-particle mode functions Xk are defined by 
(lA-lOp . Defining the power spectrum in terms of the 
two-point function in the usual manner 



5(3)(k-|-k')|fp0 (A-13) 



we can extract the power in re-scattered (p modes. 

Alternatively one could compute the power spectrum 
of re-scattered inflaton modes using the Schwinger's "in- 
in" formalism [2^ which was implemented to compute 
cosmological perturbations by Weinberg in ^30] . We have 
verified that the tree level contribution to obtained 
using this formalism reproduces our result dHS). O ur 
approach is analogous to computing the cosmological per- 
turbation from the field equations using the Seery et al. 
approach [3l|. The consistency of this method with the 
in-in approach at tree level is in accordance with the gen- 
eral theorem of f32l . 



A. 3. Renormalization 

To compute the spectrum of re-scattered i5(/)-particles 
we simply need to insert the solution (jA-8[) into (jA-12[) 
and evaluate the integrals. However, there is one sub- 
tlety. The resulting power spectrum is formally infinite, 
moreover, it contains the effect of both particle produc- 
tion as well as vacuum fluctuations of the x field. We 
are only interested in the re-scattered 5(j) which are due 
to particle production, thus, we need to subtract off the 



10 



contribution due to nonlinear 54> production by x vacuum 
fluctuations. 

To properly define the two-point function of 5(j) we need 
to renormalize the four-point function of the gaussian 
field X- ^ warm-up, let us first consider how to renor- 
malize the two point function of the gaussian field x- We 
use the following scheme 

(A-14) 

where (^^^ (^i)Cfe2 (^2))in is the contribution in the absence 
of particle production, computed by simply taking the 
solution (|A-5|) with ak — I, f3k = 0- More explicitly, for 
the case at hand, we have 



(X'(i,x)),, 



(27r)3 
(x'(i,x))-5 



\xlit)\ 



(A-15) 



where Sm is the contribution from the Coleman- Weinberg 
potential. This proves that our prescription reproduces 
the one used in \W^ . 

Having established a scheme for remormalizing the 
two-point function of the gaussian field x it is straight- 
forward to consider higher order correlation functions. 
We simply re-write the four point function as a product 



J 



of two-point functions using Wick's theorem. Then each 
Wick contraction is renormalized as above. Applying this 
prescription to (jA-12p amounts to 

X J dt'dt"t't" sm[nk^{t - t')]s\n[nk^{t - t")] 

X j d'k' [Xk.-k'it')xl-k'it")~h.-k'it')fl-k'it")] 
x[xk'it')xl,it")-h'it')n'{t")] (A-16) 
where fk{t) are the adiabatic modes defined in (|A-4p . 



A. 4. Spectrum of Re-scattered Modes 

Let us now proceed to compute analytically the renor- 
malized spectrum of re-scattered inflaton modes by 
inserting the solutions (jA-8|) and (|A-4[) into (jA-16p and 
carrying out the integrations. The computation is te- 
dious but straightforward since the time and phase space 
integrals factorize. We have relegated the technical de- 
tails to appendix B and here we simply state the final 
result: 
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fc3 



-7i-fc^/(2fcf) 

2V2 



7r„ 
-Re 
4 



3\/3 5^5 



^2inkt-ini/i2ki)-iTT /2 p 



^ [1 - cos{nkt)f 



[1 — cos(r2fet)] Im 



(A-17) 



Equation (jA-17P is the main result of this appendix. The the solutions (|A-8p and (|A-4p into (|A-16p we find the 
"form factor" F(k, t) is given explicitly in appendix B. result 



APPENDIX B: Detailed Computation of 



In this appendix we discuss in some detail the technical 
details associated with the computation of P^. Inserting 
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(Tk'nk-k'nw / dt'dt" sin [Q.k{t - t')\ sin [17fe(i - t")\ cos^ 



X j dt'dt" sin - t')] sin [f7fe(t - i")] sin^ 



{Kt'f iKt") 



ii\2-\ 



[Kt'f , (Kt"f 



d^k' 



rik-k' V^Vl + nk' + nk'y/nk-k'\/l + Uk-k' 
X I dt'dt" sm[ak{t-t')]sm[nk{t-t")]sin ^* ' > 



{k.t'f {k,t"f 



(B-1) 



We consider the time and phase space integrations sepa- 
rately. 



B.l. Time Integrals 



AU the time integrals appearing in (jB-l[) can be written 
in terms of two functions which we call /i, 12- These 
involves are defined as 

h{k,t) = / dt'sin[f7fe(<-t')]e^^''**')' (B-2) 



h{k,t) = 



[ dt'sm[nk{t-t')] (B-3) 

^0 



First consider Ii. It is useful to factorize the answer 
into the product of the stationary phase result (valid for 
ki,t > flk/{2ki,) > 1) and a "form factor" F{k,t) as 
follows: 



h{k,t) 
F{k,t) 



(1 + 6-2^"'=*) erf 



2 k^ 



(B-4) 
(B-5) 



'='erf 



The form factor F{k, t) has a complicated structure. We 
have illustrated the qualitative behaviour of this function 
in Fig. [5] taking ^Ik/ki, — 5 for illustration. 

Next, consider the characteristic integral I2, eqn. 
(|B-6[) . This integration is trivial: 



l2{k,t) = — [l-cos{nkt)] (B-6) 
life 

Now we will show that all the time integrals appearing 
in (|B-ip can be reduced to combinations of the charac- 
teristic functions Ii and 12- First, consider the first line 




FIG. 8: The behaviour of the function F{k, t) as a function 
of t. For illustration we have set Q.^ = 5fc*. 



of (|B-ip where the following integral appears: 

j dt'dt" sin [Vlkit - t')] sin [^kit - t")] 
-{k^t'f {k^t"f- 



2 2 

|/l(fc,t)|2 , /2(fc,t)2 



= ^|F(fc,t)p + ^[l-cos(f]fet)f (B-7) 



Next, consider the time integration on the second line 
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of mAl : 



dt'dt" sin [^lk(t - t')] sin [^lk{t - t")] 



terchange of a and 6. 

The phase space integral in the first line of (|B-1[1 is 
computed by a trivial application of the identity (|B-10p : 



X sm 



2 2 



d k Uk-k'nk' = 



4< 



(B-11) 



rRe 



■F{k,ty 



+ ^ [1 - cos(r!fei)]^ 



Finally, we consider the time integration on the third 
line of (foil : 

J dt'dt" sin [Qkit - t')] sin [f7fe(i - t")] (B-9) 
"(M')' , {k*t"f] r(fc.t')' 



The remaining integrals cannot be obtained exactly 
in closed form because they contain terms like VI + '^fc' 
where the gaussian factors appear under the square root. 
However, because < 1 it turns out to be a very good 
(-^"^) approximation to replace ^/l + nu' = 1 + Uk' 12. (We 
have checked numerically that the error induced is less 
than a few percent.) Let us now proceed in this manner. 
The phase space integral on the second line of (|B-ip is: 



d^k'y/nk-k'nk'\/l + Uk-k'V^ + nfe' 



X sm 



2 2 

= Im[/i(fc,i)/2(fc,<)] 



2kM 



k 



[1 — cos(i7fe<)] Im 



'F{k,t) 



d^k' 



exp 



1/2 1/2 1 3/2 1/2 1 1/2 3/2 



2V2 



exp 



(B-12) 



B.2. Phase Space Integrals 

Throughout the calculation integrals of the following 
form appears frequently: 



d k Ti^ ^/Tif^/ 

J d^fc'exp [-a7r|k - kf //fc^] exp [-67r|kf //fc^] 

(B-10) 



Finally, consider the phase space integral on the third 
line of (|fo|) : 



d^k' 
d^k' 



Uk-k'y/nJ^V^ + nk' + nk'^/uk-k'^/l + nk-k' 



1/2 , 1/2 



3/2 



3/2 



+ -nk-^k'nt^, +2*^^'^-^' 



(a + 6)3/2 



exp 



ah irk 



a + b kl 



4V2 f 7rfc2 



3V3 



exp 



ikl 



2V2 f Snk^ 
: exp 



5kl 



(B-13) 



This formula is valid when a, b are positive real num- 
bers. Notice that this expression is symmetric under in- 



Assembling the various results presented in this ap- 
pendix one arrives straightforwardly at the result (|A-17p . 
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